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ABSTRACT 

We present the logA^-log5 and spatial distributions of X-ray point sources in seven Galactic Bulge (GB) fields 
within 4° from the Galactic Center (GC). We compare the properties of 1159 X-ray point sources discovered 
in our deep (100 ks) Chandra observations of three low extinction Window fields near the GC with the X-ray 
sources in the other GB fields centered around Sgr B2, Sgr C, the Arches Cluster and Sgr A* using Chandra 
archival data. To reduce the systematic errors induced by the uncertain X-ray spectra of the sources coupled 
with field-and-distance dependent extinction, we classify the X-ray sources using quantile analysis and estimate 
their fluxes accordingly. The result indicates the GB X-ray population is highly concentrated at the center, more 
heavily than the stellar distribution models. It extends out to more than 1.4° from the GC, and the projected 
density follows an empirical radial relation inversely proportional to the offset from the GC. We also compare 
the total X-ray and infrared surface brightness using the Chandra and Spitzer observations of the regions. The 
radial distribution of the total infrared surface brightness from the 3.6 band jim images appears to resemble the 
radial distribution of the X-ray point sources better than predicted by the stellar distribution models. Assuming 
a simple power law model for the X-ray spectra, the closer to the GC the intrinsically harder the X-ray spectra 
appear, but adding an iron emission line at 6.7 keV in the model allows the spectra of the GB X-ray sources 
to be largely consistent across the region. This implies that the majority of these GB X-ray sources can be of 
the same or similar type. Their X-ray luminosity and spectral properties support the idea that the most likely 
candidate is magnetic cataclysmic variables (CVs), primarily intermediate polars (IPs). Their observed number 
density is also consistent with the majority being IPs, provided the relative CV to star density in the GB is not 
smaller than the value in the local solar neighborhood. 
Subject headings: Galaxy; bulge — X-ray; binaries — X-ray; population 



1. INTRODUCTION 

The Chandra X-ray Observatory has opened a new era in 
studies of the X-ray source population in the Galactic Bulge 
(GB). A series of shallow and deep Chandra observations in 
the Galactic Center (GC) region have re vealed ~ 1000 X' 
ray point sources in a 2° x 0.8° region (IWang etal.ll2002h 
and 2357 X-ray point sources in a 17' x 17' region around 
the Sgr A* (Muno et al. 2003, hereafter M03). An additional 
^ 2000 sources found in the Bulge Latitude Survey (BLS, 
two 0.8° X 1.5° regions) provide the initial results for th e 
latitude distribution of the GB sources (iGrindlav et al.ll2009h . 
The X-ray luminosities and relatively hard spectra ruled out 
that the majority of the GC X-ray point sources are normal 
stars, active binaries, young stellar objects, or quiescent low 
mass X-ray binaries (qLMXBs) (M03). From the lack of real 
matches between the bright i nfrared (IR, K < 15) and X-ray 
sources in the Sgr A* field, Laycock et al.l (l2005 l hereafter 
L05) concluded that high mass X-ray binaries (HMXBs) can- 
not account for more than 10% of the X-ray sources in this 
region. While the leading candidate that fits the properties 
of thes e X-ray sources is now magnetic cataclysmic variables 
(CVs) (iMuno et al.l2004 L05), the relatively hai'd X-ray spec- 
tra of some of the most recently discovered qLMXBs imply 
qLMXBs could be misrecognized as C Vs and be more com- 
mon in the GB than th ought in the past ( Wijnands et al. 20051 
iBogdanov et al.ll2b05h . Infrared (IR) searches for the counter- 
parts oftiiese_GBX;ray sources have been actively pursued 
(e.g. iMuno et aP (l2005h ). but the exact nature of the majority 
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of the sources is still elusive due to high obscuration by dust 
and source confusion by the high star density. 

We have conducted a series of deep (100 ks) Chandra ob- 
servations of three low extinction Window fields - Baade's 
Window (BW), Stanek's Window (SW; Stanek 1998) and the 
"Limiting Window" (LW) - neai- the GC (©. These Win- 
dow fields allow us to observe the GB X-ray population and 
their Galactic radial distribution with minimal obscuration by 
dust. We have discovered 1159 X-ray point sources in these 
fields. We compare their distributions with X-ray sources in 
other GB fields - the Sgr B2, Sgr C, Arches Cluster and Sgr 
A* fields. We present a new approach using quantile analysis 
(® to minimize the systematic errors in flux estimation, to 
classify sources by their X-ray spectral types and investigate 
their radial distribution. We compare the X-ray distribution 
with the known models of the stellar distribution (® and in- 
vestigate the natur e of th e X-ray population (©. See also 
Ivan den Berg et all ( l2009l) . where we put constraints on the 
nature of the X-ray source populations from the optical point 
of view, using HST observations of the Window fields taken 
simultaneously with the Chandra exposures. This work is part 
of our Chandra Multi-wavelength Plane (ChaMPlane) Survey 
designed to measure the space density and probable nature of 
the low-lum inosity accretion-powered sources in the Galaxy 
(IGrindlav et al. 2005). 

2. OBSERVATIONS AND DATA ANALYSIS 

We performed Chandra/ AClS-l observations of BW on 
2003 July 9 (Obs. ID 3780), SW on 2004 February 14/15 
(Obs. ID 4547 and 5303), and LW on 2005 August 19/22 and 
October 25 (Obs. ID 5934, 6362 and 6365). Due to techni- 
cal constraints, the SW and LW observations were segmented 
into a few pointings, which we stacked for further analysis. 
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Table [T| summarizes the observational parameters and X-ray 
source statistics of the Window and other GB fields analyzed 
in this paper For the Sgr A* field, we use the results from 
a 100 ks observation (Obs. ID 3665) for easy comparison 
with other GB fields that were observed with similar expo- 
sure times, and we have also stacked 14 observations from 
the archive, totaling 750 ks exposure. 

We have analyzed the data as a part of the ChaMPlane sur- 
vey. For uniform analysis of all the ChaMPlane fields, we 
have developed a series of X-ray process ing tools, raainly 
based on version 3.4 of the CIAO package dHong et alj|2005l 
hereafter H05)~. After initial screening of the CXC level-2 
data (e.g. select the events in good time intervals during which 
the background fluctuates < 3(T above the mean level), we de- 
tect X-ray p oint sources with a wavelet algorithm (wavdetect, 
iFreeman et a l. 2002) with a significance threshold of 10"^. 
The wavdetect routine is run on each individual observation 
and the stacked data set if available. Multiple observations 
are considered stackable (the SW, LW and Sgr A* fields here) 
if the aimpoints are on the same detector (ACIS-I or ACIS-S) 
and they are within 1' of each other 

In H05, we used source detections in the broad {Bx'- 0.3- 
8.0 keV) band. We now also incorporate source detections in 
the soft (Sx: 0.3-2.5 keV) and hard (Hx: 2.5-8.0 keV) bands 
in addition to the broad band. We establish a unique source 
list by cross matching the three detection lists based on the 
relative distance of possibly identical source pairs (the closest 
pairs) in the three different bands. The relative distance {dr) 
of two sources is defined by the ratio of the source distance 
to the quadratic sum of the positional errors. Note that there 
is no astrometric offset among the images in the three bands. 
The positional errors of sources are calculated by an empirical 
formula based on the MARX^ simulations (Eq. 5 in H05). 

Establishing a unique source list is straightforward in rela- 
tively un-crowded fields such as the Window fields, but it can 
be tricky in heavily populated fields and in very deep expo- 
sures such as the stacked Sgr A* dataset. Fig. [T] shows the 
distributions of the relative separations of nearest-neighbor 
source pairs among the three detection bands, Sx, Hx, and 
Bx- As examples, we compare the 100 ks observations of 
the Stanek Window and Sgr A* field, and the 750 ks stacked 
dataset of the Sgr A* field. A source detection in each band 
contributes two pairs to the distribution, one from each of the 
other two bands. 

The bimodal shape of the distribution indicates two types 
of the pairs contribute to the distribution: one type consists of 
the truly identical sources detected in different bands and the 
other consists of random pairs of unrelated sources. The dis- 
tribution of the random pairs can be estimated by introducing 
an arbitrary astrometric offset in the source position between 
the detection bands. The (blue) dashed-dotted line in Fig [T] 
shows such an example (1' offset in both R.A. and Dec), the 
shape of which closely resembles the right side of the bimodal 
distribution of original sources. The slight excess over the 
original distribution is due to the real pairs being transformed 
into new random pairs by the positional offset. 

After visual inspection of the raw images and the distribu- 
tions of the relative distances in Fig. [T] we use a simple cut 
(dr <2.0, red vertical line) for establishing the unique source 
list. The cut is sufficient for identifying virtually 100% of the 

- Some of the fields were processed by the tools based on version 3.1 of 
the CIAO package, but the difference between two versions is minimal. 
3 http://space.mit.edu/CXC/MARX/. 



unique sources detected in multiple bands. From the distribu- 
tions of the relative distances of the random pairs, we estimate 
that the false random matches surviving the cut ranges from 
10 to 25 for the 100 ks observations. The corresponding num- 
ber of independent sources that might have lost by the false 
random matches ranges from 5 to 10 1%) for the 100 ks 
observations and about 100 for the stacked Sgr A* field (~ 
3%). 

When multiple detections in the three bands are identified 
as a unique source, we select the one with the smallest po- 
sitional error for the unique source list. We note that the fi- 
nal source position and error could be derived from a form of 
weighted average of astrometric properties of multiple detec- 
tions. However, the Sx and Hx band detections are not entirely 
independent from the Bx band detection. Therefore, in order 
to avoid unnecessary complication in the analysis, we simply 
take the astrometric (and photometric) properties of the source 
with the smallest positional errors, which is the detection with 
the highest significance among the three bands. 

As a sanity check, we compare our detection with the avail- 
able catalogues in the literature. M03 provided a catalogue of 
the 2357 X-ray point sources discovered i n the 690 ks expo- 
sure (626 ks of GTIs) of the Sgr A* field. lMuno etaH (l2006l 
hereafter M06) for the 397 point sources i n the 100 ks expo- 
sure of the Sgr B2 field (Obs. ID 944) and lWang eta"n(l2006l 
hereafter W06) for the 244 X-ray point sources in the 100 ks 
exposure of the Arches Cluster (Obs. ID 4500). The majority 
(~ 85-90%) of these sources are also detected in our anal- 
ysis or vice versa if our source list is shorter than theirs. A 
small fraction of the sources are missing due to many subtle 
differences in the detection methods and the selection criteria 
for the lists such as the detection energy bands (0.5-8, 0.5- 
1.5 and 4-8 keV in M03 or 1^, 4-9 and 1-9 keV in W06), 
the pointing or GTI selections (e.g. Obs. ID 1561 was not 
included in our analysis in the stacked Sgr A* field). In ad- 
dition, in the case of W06, their high significance threshold 
(10~^ except for the inner 2' x 2' region, about 10 counts in 
their broad band, 1-9 keV) is one of the main reasons for 
the difference (244 vs. 423) in the total number of the source 
detections. M03 and M06 list all the source detections includ- 
ing ones with negative counts in their broad band (0.5-8 ke V). 
Our source list includes the detections with net counts > 1 in 
the Bx band (0.3-8 keV). 

After source detection, we perform aperture photometry on 
each source to extract the basic source properties such as net 
count and net count rate in the conventional energy bands (Sc- 
0.5-2.0, He: 2.0-8.0 and Be'. 0.5-8.0 keV) and energy quan- 
tiles in the broad band (Bx'- 0.3-8.0 keV). For the sources that 
fall near other sources, we carefully revise the aperture of the 
source regions by excluding overlapping sections to minimize 
the contamination from the neighbors (H05). Table [2] lists a 
part of the source catalog with selected source properties used 
in this paper. The complete list for the Window and other GB 
fields is available in the electronic edition. 

3. FLUX ESTIMATION BY QUANTILE CLASSIFICATION 

In order to compare source distribution in various regions 
of the sky with diverse extinctions, it is necessary to correct 
for the interstellar absorption and use the unabsorbed source 
flux of individual sources. However, such a calculation is 
not trivial for X-ray sources with diverse spectral types as 
found in Galactic plane fields since faint sources are unsuit- 
able for spectral fitting. Moreover, the relatively large ex- 
tinction in the GB fields and its usually-unknown field-and- 
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distance dependent variation make it difficult to identify the 
underlying X-ray spectral model (e.g. power law vs. thermal 
Bremsstrahlung, etc). An inaccurate assumption of the spec- 
tral model when estimating flux introduces systematic errors 
that often exceed the stati stical errors. We therefore employ 
quantile analysis (iHong et al. 2004 . hereafter H04), which is 
relatively free of the count-dependent bias inherent in X-ray 
hardness ratio or X-ray color analysis, and so provides a bet- 
ter measure for classifying X-ray sources in the GB fields. In 
the following, the energy quantile /i j corresponds to the en- 
ergy below which x% of the counts are detected, e.g. £50 is 
the median energy. 

3.1. Quantile Analysis 

Fig. |2] shows the quantile diagrams of the X-ray sources 
> 3 in Bx) in the selected GB fields overlaid with grids 
for a simple power law model (PL, solid lines) and a power 
law plus an iron emission line (PL-nFe, dashed) at 6.7 keV 
with 0.4 keV equivalent width (EW, see ^3.2l for the motiva- 
tion of the line choice). The lower right panel also includes a 
grid for the thermal Bremsstrahlung model (TB, dotted). The 
here is calculated based o n the statistica l error s {uc) using 
small-num ber statistics from Gehrels et alJ (1986) (see also 
iKim et an (i2004)). The difference in the model grids between 
PL and PL-nFe is only evident in the highly absorbed or spec- 
trally very hard section of the diagram (the right side) because 
a small iron line (> 6 keV with < 1 keV EW) does not make 
a noticeable difference in three quantiles of the soft sources. 

Relatively insensitive to the extinction, the sources around 
(x,y) = (-0.9, 1.6) are present in every field and they appear 
unabsorbed and intrinsically soft regardless of the assumed 
model class (PL, PL-nFe or TB). Foreground thermal sources 
such as coronally active stars fit the description. The loca- 
tion of relatively hard sources in the diagram varies with the 
field extinction. In the Window fields, the hard sources are 
relatively unabsorbed, but on approach to the GC, there is 
an increasing trend in the source number with both the aver- 
age absorption and the intrinsic hardness, when compared to 
a simple power law model. For instance, in BW most sources 
have power law photon index (F) > 1 and Nmi < 1, whereas 
in LW many sources lie in F < 1 and Nmi > 1 . In the Sgr 
A* field and the rest, most of the hard sources are heavily ah- 
sorbed with F < 1 and Nmi 1 on average, well separately 
from the foreground sources. We explore this more in Section 
3.2. 

The quantile diagrams nicely illustrate the spectral diversity 
of the X-ray sources in the GB fields, but poor photon statis- 
tics also contributes to the scatter To reduce systematic errors 
caused by poor statistics in assigning spectral types while al- 
lowing the spectral diversity of the sources in each field, we 
divide the diagram into three groups as shown by the (grey) 
lines originating at (-0.5, 1.3). The left section represents 
most foreground thermal sources (Gl: soft group), the mid- 
dle section most unabsorbed accreting sources (G2: medium 
group), and the right section the absorbed thermal or accret- 
ing sources (G3: hard and absorbed group). The division be- 
tween Gl and G3 is devised to be somewhat robust'* against 
variations in detector response between ACIS-I and ACIS-S 
(see H04; H05); or induced by gradual loss of low energy re- 
sponse. The final results (e.g. logA^-logS' distributions) are not 
sensitive to small changes of the group boundaries (e.g. shift- 

* Under the PL model, the boundary of the Gl and G2 groups stays in 
roughly between r=2 and 3. 



ing the boundaries by ~ 0. 1 in x or y). The mean quantiles for 
each group (marked by 0) are calculated by the stacked pho- 
tons of the sources in the group with S/N >3 and net counts 
< 1000 (to avoid being dominated by a few bright sources) in 
Bx- For a given model class (e.g. PL), we estimate the spec- 
tral model parameters (e.g. F and A^h) of the sources in each 
group using the mean quantiles. 

3.2. Spectral hardening vs. radial offset from GC 

Table [3] summarizes the group mean quantiles of the G2 
and G3 sources and corresponding model parameters under 
the PL and PL-nFe models. The G2 sources in the high ex- 
tinction fields are omitted in the Table since they are mostly 
foreground sources. For comparison, the table also shows the 
model parameters estimated from spectral model fits. In order 
to increase photon statistics, we stacked the spectra of sources 
within a group with net counts < 1000 for spectral fitting^ and 
the spectra were binned to have at least 40 counts in each bin. 
The model used is a power law plus an iron emission line 
for which we have chosen the 6.7 keV Fe XXV He- a line 
because it has also been observed in the spectra of the X-ray 
point sources in the deep survey of the Sgr A* field (M03), the 
shallow survey of the GC strip (Wang et al. 2002), and other 
parts of the Galactic plane (iEbisawa et al.ii2005i) . The 6.4 keV 
neutral iron line is also present in some sources of the GB 
fields, but it is generally m ore prominent as unresolved dif- 
fuse emission (IWang et al.ii200 2). Note our aperture photom- 
etry is designed to minimize possible contamination of the 
diffuse emission through background subtraction (H05). We 
have chosen the 0.4 keV EW for the line in quantile analysis 
(and the logA^-logS' distribution later) because it lies in the EW 
range estirnated by spectral model fits on the G3 sources and 
iMuno et aP (l2004l) found a similar value (~ 0.4 keV) for the 
bright sources in the Sgr A* field. 

Under the PL model, both groups show a trend of increasing 
hardness of the intrinsic spectra on approach to GC (i.e. F > 1 
for the Window fields vs. F < 1 for other GB fields). This 
apparent trend can be attributed to a few factors. 

In the case of the G2 sources in the Window fields, the 
group contains a large number of foreground coronal sources 
and background AGN in addition to the GB sources. For in- 
stance, in the BW and SW, about 70% of the sources with 
SjN > 3 in the He band are background AGN (see Fig. |5]in 
Therefore, the group is perhaps too contaminated for the 
apparent trend to be taken for real. 

The similar trend in the G3 group appears to be more re- 
alistic. However, comparing the PL-nFe model grid with the 
PL model grid in the quantile diagram suggests this trend can 
be an artifact of using the PL model at least in part. Indeed 
the trend is alleviated under the PL-nFe model as shown by 
the model parameters estimated by both quantile analysis and 
the spectral fits (TableO. Poor statistics in the G3 sources of 
BW and SW does not allow any meaningful constraint of the 
iron line emission in the spectra, but the stacked spectrum of 
the G3 sources in the L W does show a clear hin t of the 6.7 
keV line (Fig. [3] see also ' Revnivtsev et"al ( l2009l) ). The rela- 
tively weak line feature (EW ~ 0.17 keV) in the LW can be 
explained by the relatively large contribution of AGN in the 
group compared to other GB fields (see Fig.|5]in ®. 

^ One can use a spectral model fit on individual sources with net counts > 
200-300, but in order to establish more reliable statistics for the presence of 
the line emission for the group, we also stack moderately bright sources with 
net counts up to 1000. 
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With the inclusion of a 6.7 keV Une, the power law index 
(F) becomes largely consistent across the GB fiel ds w ith the 
possible exception of the Sgr B2 or C field (see 35.4b . The 
result indicates that the galactic X-ray sources in the LW field 
may be the same type of sources as seen in the other GB fields 
closer to the GC. If true, the GB X-ray sources indeed extend 
out to at least 1.4° from the GC. 

3.3. Flux Estimates 

Based on the group model parameters, we estimate the 
source flux in the conventional energy bands, using the in- 
strument response files at the aim point and scaling them 
for source position on the detector by the exposure map 
(H05)^. For the stacked data, we use the average response 
files weighted by the exposure of each observation. 

The quantile diagram can assign the model parameters (e.g. 
F=1.7 vs. 1.0) appropriate to a given model class for the 
sources, but it cannot determine which model class (e.g. PL 
vs. TB) is right for the sources. A certain model can only be 
ruled out when the derived values of the parameters are un- 
physical or with external information (e.g. optical identifica- 
tions). In order to estimate the systematic errors arising from 
the improper choice of the model class, we compare flux esti- 
mates under three different model classes: PL, PL-nFe and TB. 
In order to see the significance of the difference among these 
models, Fig.|4]compares the conversion factor of count rate to 
unabsorbed flux for sources near the aim point in each group 
under the three model classes. In the case of the He band, the 
difference between the model classes is very small, but in the 
Sc band, the conversion factor can differ by more than a fac- 
tor of 10. We take the largest difference in the flux estimates 
among the three model classes as the systematic errors (ct,) 
and compare them with the pure count-based statistical errors 
(<Tc). For the flux estimates in the He band, (t, ^ 20-30% and 
we get (Tj < CTc for ^ 87% of the sources with S/N > 3, and 
cTj < (Tc for ~ 62% even with S/N > 5. In the Sc band, CTj can 
be larger than 1000%, and we get cr, < dc only for 28% with 
S/N > 3 and for 16% with S/N > 5. 

This exercise does not explore all the possible model 
classes, but the results indicate the He band flux estimates in 
this method are robust and relatively insensitive to the choice 
of the model classes. However, the Sc band flux estimates can 
be dominated by the systematic errors arising from improper 
selection of the spectral model. The fundamental difference 
between the Sc and He band is that the Sc band is very sen- 
sitive to the range of interstellar absorption in the GB fields, 
~ 1021-23 cni"^, while the He band is not. In the following, we 
limit our discussion to the He band results using the sources 
with S/N > 3 in He ('x's and triangles in Fig. |2ll, which are 
likely to be GB sour ces (and AG N) rather than the foreground 
sources. See also v an den Berg et al. (2009 ) for the spectral 
choices for the flux estimates in Bx- If the number of counts in 
the stacked spectrum of a quantile group (Gl, G2 or G3) were 
large enough, a spectral fit would be better suited for deter- 
mining the underlying spectral model and its parameters, but 
a fit can also leave ambiguity over the correct spectral model. 
Since in the He band the difference driven by the model class 
is less significant than the statistical errors, we simply use the 
model parameters estimated by quantile analysis. 

4. SOURCE DISTRIBUTION 



4. 1 . Eddington and Malmquist Biases 

In order to explore the effect of the Eddington bias (EB), 
which makes the faintest sources appear brighter than they re- 
ally are, we simulate three spectral types of sources (one for 
each group) based on the group mean quantiles for each field. 
Using the MARX simulation code, we generate the sources 
with net counts {Bx) from 5 to 400, using a distribution 
to cover the wide count range efficiently. We scatter 200 - 
250 of these sources randomly over the real events and apply 
the regular analysis procedure. We repeat the procedure 1000 
times. The fake sources are not allowed to overlap each other 
but they can fall on top of the real sources. The results indi- 
cate that the EB is noticeable in the sources with < 10 counts, 
which can appear as bright as 15-20 count sources, depend- 
ing on the field (or up to ^ 30 count sources for the stacked 
Sgr A* field). Since we consider sources with S/N > 3 in the 
He band, which corresponds to > 16 - 20 counts at least (> 
30 counts for the stacked Sgr A* field), we expect EB is not a 
major contributor to the errors of the following distributions. 

The Malmquist bias (MB) is due to the exposure dependent 
volume (depth) coverage. The MB is usually a concern for 
luminosity distributions but not for logA^-log5 distributions in 
the apparent (detected) flux space. However, the logA^-logS' 
distributions in the unabsorbed flux space can be subject to 
the MB when strong interstellar absorption limits the depth of 
the view, underestimating the true distribution. Therefore, the 
faint end of the logA^-log5 distribution can be lower than the 
true distribution and so the MB counteracts the EB to some 
extent. 

With > 100 ks exposure, all the sources with an unabsorbed 
flux > lO"''* ergs cm~^ s"' can be detected at the far side of the 
Galaxy with > 3 in He under the assumption of the total 
integrated absorption^ of A^h 12 x 10^^ cm~^. Therefore the 
MB is not a concern for sources with > 10"'^ ergs cm"^ s"' 
(or > 5 X 10~'^ ergs cm"^ s"' for the Window fields) under 
the assumption of a power law wth r=1.0 for the X-ray spec- 
trum. This does not mean we can access X-ray sources of a 
certain luminosity uniformly all the way through the Galaxy. 
For instance, the unabsorbed flux of 5 > 10"''* ergs cm"^ s~' 
coiTesponds to 8 x 10"" ergs s"' at the GC (8 kpc, A^h 
- 6 X 10^2 cm-2) and Lx>l x lO^^ ergs s"' at 20 kpc (A^h 
~ 12x10^^ cm~^). The situation is a bit more tricky since 
under the quantile group method we assign fixed spectral pa- 
rameters with a fixed A^h value for the X-ray spectra of all the 
sources in each group, which in fact have a diverse A^h distri- 
bution (e.g. the G3 group in the Sgr A* field). However, the 
sources with an unabsorbed flux of 5 > 10"''* ergs cm"^ s"' 
in this method are free of the MB for > 100 ks exposure. 
The MB can be a problem for the Gl and G2 sources in the 
high extinction fields, but their contribution in the logA^-logS' 
distribution of the He band is negligible compared to the G3 
sources. 

4.2. Sky Coverage 

For the logA^-log5 distribution, we need to know the sky 
coverage as a function of flux. In order to minimize the sys- 
tematic errors associated with spectral type assignment, we 
calculate the sky coverage of each source based on the de- 
tected photon counts in the three energy bands as follows 
(Cappelluti et al. 2005; H05). For each observation, we gen- 
erate the background-only images by removing the counts in 



' The latest CIAO tools (ver 3.4 or higher) can calculate the response files 
appropriate for each source location. 



^ A ssuming the absorption to the GC to be 6 X 10^^ cm ^ jBaganoff et all 
120031) and the symmetry with respect to the GC. 
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the source regions in the image and filling the region with 
the counts using the statistics in the surrounding regions {dm- 
filtlif. At every pixel in the background images, we calcu- 
late the minimum source counts required for detection with 
S/N > 3. For the sky coverage of a given source, we take 
the sky area where the minimum counts are less than the net 
counts of the source in the band. These sky coverage values 
agree well with those expected from the simulated sources of 
three spectral types using the MARX 04.1l i. On average, they 
are within 10% for the cases with S/N > 3, which indicates 
this method accounts for the completeness as well. 

4.3. The logN-logS and Radial Distributions 

The logA^-logS' and radial distributions of the X-ray sources 
in the GB fields are shown in Fig. [5^ & b respectively. 
The logA^-log5 distribution was computed using sources with 
S/N > 3m the He band, under both the PL and PL-nFe mod- 
els described in ^ with the latter result plotted in Fig. |5^. 
The source number-density values plotted against angular dis- 
tance from Sgr A* in Fig. ^ are projected from the logA^- 
log5 distributions at the flux value indicated by the vertical 
grey line (S > = 1.5 x 10"'"^ ergs cm~^ s"') in Fig. |5^. As 
seen in Fig.jSj), the total source densities under the PL model 
are slightly lower than the same under the PL-nFe model in 
the high extinction fields, while both distributions under these 
models are nearly identical in the three low extinction Win- 
dow fields. 

For clarity, we define the statistically robust section of each 
distribution in Fig.|5^ and emphasize it with a solid line. This 
"solid section" is defined to contain contributions from at least 
10 or more sources, which set the upper limit of the range (e.g. 
So ^ lO"'-' ergs cm"^ s"' for Sgr C). The lower limit is set by 
the flux value at which the sky coverage of the contributing 
G3 sources is greater than 50% of the maximum sky coverage 
i.e. the full field of view (FoV) (e.g. So ^ 10"'"* ergs cm"^ s"' 
for Sgr C). In this way, we avoid the statistical bias or fluctua- 
tion due to either low source statistics at the bright end of the 
accessible flux range or limited sky coverage at the faint end 
of the range. The portion of each distribution not meeting the 
above criteria is dotted. 

The slope (a) of the logA^-log5 distribution is calculated by 
a power law fit (A^ oc to the solid line section of the dis- 
tribution. The y-axis error of the distribution is given by the 
quadratic sum of the statistical error (shown in the figure) and 
a constant systematic error ('^ 20%, the difference between 
the PL and PL-nFe model). As expected for the narrow FoV 
of ACIS-I observations, the slopes of the logA^-logS' distribu- 
tions are largely consistent with the -1.5 slope within ~ 2a 
except for the stacked Sgr A* field, which shows a hint of 
the actual deviation 6<t) from the -1.5 slope. Note the cal- 
culated slopes are only for guiding purpose, and they should 
not be taken seriously for representing the population since a 
simple power law does not fit some of the distributions very 
well. 

The AGN disti'ibution is taken from ' Kim et al.l (1200 7^. us- 
ing a power law model with F = 1.7 for the X-ray spectra. 
The (orange) dashed line indicates the reduced AGN popula- 
tion that can be seen through the high extinction fields such as 
the Sgr A* field, since the unabsorbed flux is corrected for the 
average absorption of the X-ray sources, mostly Galactic and 
centered around the GC, which should be about half of the to- 
tal absorption for the AGN. For simplicity, we correct another 

^ http://cxc.harvard.edu/ciao/ahelp/dmfilth.htmI 



Nu = 6x 10^^ cm"^ for the AGN seen in the high extinction 
fields. 

For the Sgr A* field, we plot the results from both the 
stacked data (black) and the 100 ks exposure (grey) in Fig. 5a 
and use only the stacked data in Fig. 5b. The spectral models 
from the 100 ks exposure are used for both data sets for fair 
comparison with other fields and to avoid any spectral param- 
eter driven variations between two exposures for the Sgr A* 
field. The distribution of the stacked data is ~ 40% higher at 

= 10"'"^ ergs cm"- s~' than the same for the 100 ks expo- 
sure. A few factors such as the MB^ are responsible for the 
difference, but the main cause of the difference is suspected 
to be the X-ray variability of the sources. The stacked data 
(750 ks) simply have a better chance of detecting the sources 
or catching high flux states of the sources than for the shorter 
exposure (100 ks). For instance, the 20 brightest sources in 
the He band in the 100 ks observation of the Sgr A* field are 
found to be about 30% brighter on average in the stacked data 
set, and five of the 20 brightest sources in the stacked data 
were not detected in the 100 ks observation. This variation 
qualitatively agrees with the change seen in the logA^-logS' dis- 
tributions, but the diverse nature of the X-ray variability and 
duty cycles makes it hard to quantify the resulting difference 
in the logA^-logS' distributions. 

The radial distribution is generated from the sources with 
5 > = 1 .5 X lO"''* ergs cm"^ s"' , plotted over the stellar dis- 
tribution (green) and the 1/9 distribution (blue dashed). Both 
the stellar and 1/6 distributions are averaged over the ACIS-I 
FoV (17' X 17') of the GB fields. The So value for the radial 
distribution is chosen as a compromise between having suf- 
ficient source statistics in the Window fields (5o < 2 x 10"'"* 
ergs cm"^ s"') and avoiding statistical biases in the high ex- 
tinction fields {So ^ 10~^* ergs cm"^ s"'). The curve resulting 
under the PL-i-Fe model (black solid) is more centrally con- 
centrated around the GC than the PL model (red dotted). 

The stellar distribution is derived from the Galactic stellar 
models compiled by M06, and it consists of a central spher- 
ical cluster (pa, Eq. 1 in M06), a central disk (pg, Eq. 2 in 
M06), a triaxial ellipsoidal GB (pc, Eq. 3 in M06) and a 
Galactic disk (po, Eq. 6 in M06; for the origin of the for- 
mulae, see also Launhardt et al. 2002, hereafter L02; Kent et 
al. 1991, hereafter K91). For the first three components {p^, 
Pb & Pc), we use the formula and the parameter values in 
L02 and M06. For the Galactic disk component (po), M06 
use a simple exponential form in K91 and employ 10" M© 
for the total Galactic disk mass for the overall normalization, 
but since the first three components are mainly for the stellar 
mass, we believe this is an overestimate. Therefore, we use 
a normalization that matches the l ocal stellar mass density of 
0.044 Mq pc"^ (iRobin et alJl2003h . This gives 1.8 x 10'° Mq 
for the w hole disk, which is roughly consistent with the es- 
timate by iRobin et al.1 (l2003h (2.2 x 10'" Mq)'*'. Since we 
expect both the X-ray and stellar sources are centrally con- 
centrated around the GC, in Fig. 5b we further assume that all 
detected hard Galactic X-ray sources are at a distance of 6 - 
10 kpc, which is justified given that the stellar models predict 
that > 80% of sources along the line of the sight of the GB 
fields lie in the same distance range. Note the central concen- 
tration also makes our normalization change of the Galactic 

' Note that the unstacked Sgr A* field (Obs. ID 3665) ha.s the shortest GTI 
(88 ks) among the seven fields. 

The small difference is mainly due to the difference in the assumption 
of the distance to the GC: 8 kpc for the model used here, and 8.5 kpc for 
IRobin et all <200l . 



6 



disk component less important in the outcome, but we find 
that the change makes this Galactic stellar mass model con- 
sistent with other Galactic stellar number density models (see 
95. H and Table|4|i. These stellar model components have about 
a factor of two uncertainty (M06, L02). 

The normalization of the stellar and 1/6' distributions is set 
by a simple fit to the radial-distribution curve under the 
PL-nFe model (Fig. 5b). We use the stacked result for the Sgr 
A* field. The radial distribution shows the GB X-ray sources 
are highly concentrated at the GC, more than the stellar distri- 
bution. It also shows the hard GB X-ray sources extend out to 
> 1.4° from the GC, roughly following an empirical relation 
of 1 /6' with some excess in the Arches Cluster, Sgr C and Sgr 
B2 fields. 

The excess of the X-ray source to stellar distribution near 
the GC does not appear as prominent if we use the 100 ks ex- 
posure of the Sgr A* field at S[) = 1.5 x lO"''* ergs cm"^ s"^. 
However, the trend of the relative excess of the X-ray sources 
toward the GC is present from the Sgr B2 to Sgr A* field (e.g. 
the Sgr B2 field has a deficit under the current relative normal- 
ization in Fig. 5b), and the logA^-log5 distributions of the 100 
ks exposure and the stacked Sgr A* field become more consis- 
tent at > 2 X 10"'"* ergs cm"^ s"'. Therefore, the excess of 
X-ray sources toward the GC with respect to the considered 
stellar model appears real. We explore this excess in more 
detail in Section 5.1. 

5. DISCUSSION 

We find that the number-density of the hard X-ray sources 
in the GB is significantly elevated above the AGN density out 
to at least the LW at 1.4° separation from Sgr A* (Fig.|5^). 
Furthermore, although empirical (see Fig. |6] for the compo- 
sition of the stellar population), the radial distribution of the 
hard X-ray sources roughly follows a 1 /6' relation out to this 
field (Fig. This discovery suggests that all such sources 
observed within at least ^ 200 pc of the GC belong to the 
same centrally concentrated population. The similarity of the 
stacked spectra of the hard X-ray sources in all fields from 
Sgr A* to LW, and in particular the presence of a 6.7 keV iron 
emission line, strengthens our conclusion that a single under- 
lying class of sources makes up this population. 

5.1. X-ray Source Density vs. CV density 

The current leading candidate to explain the X-ray sources 
within 20 pc around the GC is magnetic CVs or intermedi- 
ate polars (IPs) in particular (M04, L05). Recent popula- 
tion synthesis models by iRuiter et al.l (i2006t hereafter R06) 
show IPs can constitute the majority of these X-ray sources 
under the assumption that IPs span a luminosity range of 
- 3 X 10^''-5 X 10^^ ergs s"' and that they make up --2-8% 
of all CVs (see also M06 for a review of the population syn- 
thesis models for the X-ray sources in the GB). Now our ra- 
dial distribution indicates this source population extends out 
to ~ 200 pc and the hard X-ray spectra with the iron emission 
line supports the idea that IPs are the major component of the 
population. 

In this section, we compare the observed X-ray source den- 
sity with the stellar (mass) density to see if CVs, especially IPs 
can explain the majority of the detected X-ray sources. Table 
|4] summarizes the number density of the X-ray sources with 
S > 1.5 X 10"'"^ ergs cm"^ s~' in the He band and compares 
them with the stellar (mass) density. For the GB X-ray source 
densities, we subtract the expected number of AGN, which is 
145 deg~^ in the Window fields and 107 deg"^ in the high ex- 



tinction fields, from the surface density (see ^4.3l l dKim et al.l 
l2007h . 

Table |4] quotes the average stellar density over the volume 
defined by the distance between 6 and 10 kpc in the 17' x 17' 
FoV using two stellar models. Stellar model A is the same 
stellar mass model used in Fig. |5](M06, L02, K91). Since 
model A provides the stellar mass density, we use a local value 
of 0.144 stars pc~^ and 0.044 M© pc"^ to convert the mass 
densi t y to the number de nsity or vice versa (iPicaud & RobinI 
2004; Robi n et al.1l2003h . As a consistency check, we also 
comp ute the stellar density usi ng another stellar model (model 
B) by iPicaud & RobinI (|2004 . which consists of a Galactic 
disk and a Galactic Bulge. This model describes the stellar 
number density in the outer Galactic Bulge and the Galactic 
disk. So it is properly normalized at the local Galaxy (0.144 
star pc"^), but due to the lack of the Galactic nucleus compo- 
nents (pa&pb), the stellar density for the GB fields within 1° 
of the GC is underestimated. The two models agree within 
30% for the Window fields where Galactic nucleus compo- 
nents are relatively unimportant. In the following, we use 
model A for comparing the X-ray source density with the stel- 
lar density. 

The relative X-ray source to stellar mass densities in the 
seven GB fields are 0.3- 1.8 x 10~^ X-ray sources Mq at 
S > 1.5 X 10"'"^ ergs cm"^ s"' (1.1 x lO-'^ ergs s"' for sources 
at the GC, 8 kpc). The large variation of the relative density 
among the seven fields reflects the mismatch between the X- 
ray and stellar distributions - the X-ray sources are more cen- 
trally concentrated than the stellar sources. The relative X-ray 
source to stellar number densities are 0.9-5.7 x 10"^ X-ray 
sources star"' at 5 > 1.5 x lO"'"* ergs cm"^ s"', depending on 
the field. 

Now we assume IPs are 5% of all CVs (e.g. ~ 2-8% for 
the models in R06) and about 12% of IPs have the X-ray lu- 
minosity abo ve 10^^ ergs s"' (e.g. ~ 10-16% in R06, see also 
iHeinke et all (l2008h ). Then the required CV to stellar den- 
sity to explain the hard X-ray GB sources ranges from 1.6 to 
9.5 x 10"^ depending on the fields. If we assume a local star 
density of 0. 144 pc"^, these correspond to the equivalent local 
CV density of 0.2- 1.4 x 10"^ pc"-*. Considering the current 
local CV density estima t es (1 - 3 x 10"^ pc"-^) in the literature 
(see e .g. lAk et al.ll2008t iGrindlav et alll2005l; iPretorius et 
l2007h . this result indicates IPs can be the major component of 
the observed X-ray sources as long as the relative CV to stel- 
lar density in the GB is comparable to the value in the local 
solar neighborhood. 

Note there are a few caveats in this analysis. First, the radial 
distribution of the X-ray sources does not match well with the 
stellar distribution as shown in Fig. 5b. As mentioned, this is 
the reason for the large variation in the estimates of the rela- 
tive CV density. It means the stellar model we use may not 
be appropriate for scaling the observed X-ray population di- 
rectly. A solution could be found in some of the assumptions 
we have made. For instance, the (apparent) fraction of IPs in 
all CVs may not be constant across the fields. 

Second, there are large uncertainties in the model param- 
eters and various assumptions such as the ratio of IPs to all 
CVs and the fractional IPs with the X-r ay luminosity > 10^^ 
ergs s~' . For instance, according to Ritt er & Kolbl (l2003l) . the 
ratio of the known IPs to all known CVs are about 10%, but 
this is also subject to a large uncertainty due to selection bi- 
ases. Similarly there is no firm estimate of the X-ray IP lumi- 
nosity distribution to set the accurate limit for the fractional 
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IPs with the X-ray luminosity > 10 ergs s . 

Third, as illustrated in the logA^-log5 distributions of the 
100 ks and stacked data of the Sgr A* field, the X-ray vari- 
ability can change the apparent source distribution. In order 
to understand the true distribution, it is necessary to monitor 
the GB fields continuously and extract the source distribution 
from a longer exposure. Considering the X-ray variability of 
the sources observed in the Sgr A* field, the true distribution 
of the GB X-ray population in the other GB fields can be > 
20-30% higher than what has been observed in the 100 ks 
exposures. 

A more recent study by ISchodel et alJ (l2007h of the stellar 
population and mass content in the Galactic nucleus shows 
evidence for a higher mass contained in the central ~ 1 pc 
than predicted by studies and they speculate the excess may 
be due to a stellar remnant population such as black holes. 
Fig. |6] shows the composition of the stellar components {p^, 
Pb, Pc & Pd) in comparison with the observed projected X-ray 
source density. Fig.|6h shows the fit results using the stacked 
Sgr A* data (the same as Fig. [S}?) and Fig. |6j3 for the 100 
ks Sgr A* data. In order to allow the possible excess of the 
X-ray population in the Galactic nucelus relative to the ob- 
servable stellar population, we also fit the X-ray distribution 
by freeing the relative normalization parameter of the central 
spherical cluster component (pa) with respect to the rest of 
the components (c & d in Fig. |6]l. The resulting fits are sub- 
stantially better, but they require a large excess of the central 
stellar cluster component relative to the fixed normalization 
given in Section 4.3: 8.4 x more for the stacked Sgr A* data 
and 4.5 x more for the 100 ks observation. This excess is still 
larger than the mass excess 2x) of t he nucleus above th e 
observable stellar cluster estimated by Schodel et al.l (l2007l) . 
although the latter has a large uncertainty (~ 30-50%). 

5.2. Total Surface Brightness 

Since the IR surface brightness is a good indicator of the 
stellar population, a direct comparison of the IR and X-ray 
surface brightness provides an independent clue on the origin 
of the X-ray emission in the fields without influence of the 
uncertainties in the stellar model. In Table |5] we list the total 
X-ray and IR surface brightness of the seven fields. In Fig.|7] 
we compare the total IR surface brightness with the X-ray 
surface brightness after removing the contribution from the 
cosmic X-ray background (CXB). The CXB is expected to 
be (1 .6 - 1 .7 ± 0.2) x 10"" ergs cm'^ s"' deg-^ for the Win- 
dows fields and (1.3-1.4±0.2) x 10~" ergs cm"^ s"' deg"^ 
for the rest dHickox & Markevitch l '2006). For easy compari- 
son. Fig. |2] also overlays the distributions of the X-ray source 
density and what the stellar model predicts without the AGN 
component in the dotted lines. 

The X-ray surface brightness is calculated from the inte- 
grated counts in the He band, using the quantile averaged 
model spectra of the G2 group for the Windows fields and 
the G3 group for the rest (Table |3]l. In order to subtract the in- 
strumental background, we use the CXC background database 
constructed from stowed observations. For the relative nor- 
malization of the background subtraction, we use the the inte- 
grated counts in the 10.5 - 12.5 keV. The error bars in Fig. |7] 
are the quadratic sum of the statistical and systematic errors 
and the errors of the CXB estimates. For the statistical errors, 
we use the difference in the flux estimates between the PL and 
PL-nFe models. 

The IR surface brightness is calculated from the 3.6 pm 
mosaic images (1.2"pixel) generated by Galactic Legacy In- 



frared Mid-Plane Survey Ex traordinaire (GLIMPSE ) from the 
SpitzerflRAC observations (IChurchwell et al.ll2009l) . The er- 
ror is based on the statistical variation of the surface bright- 
ness in the four Chandra ACIS-I FoVs of each field. In or- 
der to estimate the unabsorbed surface brightness consistently 
with the X-ray flux estimates, we assume the A^h estimates by 
the quantile analysis in Table [3] (G2 for the Windows fields 
and G3 for the rest). We use the A^h estimates by both models 
(line A for PL and line B for PL-nFe in Fig. |7]i to explore the 
dependence of the A^h estimates. For a given A^h estimate, 
we use A^H = 1.7 9 x lO^'Ay, Ak = O.lU Ay = 0.95A^,, and 
Ai.6f,m = 0.5 Aks (iNishiyama et al.ll2009l) . 

The radial distribution of the CXB-subtracted X-ray surface 
brightness roughly matches with the distribution of the un- 
absorbed IR surface brightness under the given uncertainties 
with the possible exception of the Arches field. This indirectly 
supports the idea that they share a common origin. In partic- 
ular, the distribution of the unabsorbed IR surface brightness 
matches with the X-ray source distribution, showing higher 
concentration at the GC than the stellar distribution model. 
The origin of the discrepancy between the stellar distribution 
model and the IR surface brightness is unclear and requires 
further investigations, but note that the unabsorbed IR bright- 
ness in the high extinction fields depends sensitively on the 
extinction estimates, resulting in relatively large uncertainties. 
This dependence is also evident in the large variation among 
the ACIS-I FoVs in each field, which is at least in part due 
to the variation in the extinction across the FoVs. The radial 
distributions of the total X-ray surface brightness (Fig.|7]i and 
the X-ray source density (Fig.|5]) also shows a slightly differ- 
ent trend in a few fields such as the Arches field. For instance, 
the difference in the Arches field is due to the large contribu- 
tion of a few bright sources in the total flux. We will address 
the detail of the unresolved X-ray emission wit h respect to 
the Identified point sources in a following paper dHong et alj 
12009^ . 

5.3. Comparison with other results 

According to Eq. 5 in M03, the X-ray source density in 
the Sgr A* field is 0.60 ± 0.04 X-ray sources arcmin"^ at 
S > 1.5 X 10"'"^ ergs cm"^ s"' or 1.25 x 10"^ ph cm"^ s"' in 
the 2-8 keV range under their assumption of a power law 
spectrum with T = 0.5 and A^h = 6 x 10^^ cm"^ This is 
roughly consistent with our results, 0.86 ± 0.17 arcmin"^ at 
S > 1.5 X 10"'"* ergs cm"^ s"' from the stacked results under 
the PL-nFe model. The error is derived from the quadratic 
sum of the statistical error and ^ 20% systematic errors (the 
difference between the PL and PL-nFe model). 

Table|6]summarizes a few estimates of the specific luminos- 
ity of the Galactic X-ray point sources in the Chandra/AClS 
energy range in the literature. The range of the specific lu- 
minosity in our study is the variation among the seven fields 
under the assumption of the PL-nFe model for the X-ray spec- 
tra. Using Eq. 7 in M06 and assuming the a values in Fig. 5a 
and the number density of the X-ray sources (0.3-1.8 x 10"^ 
X-ray sources Mq at > 1 . 1 x lO-'^ ergs s"' ) in Table 4, we get 
0.5-2.8 X lO^*" ergs s"' in 2-8 keV for the luminosity 
range of 5 x lO-'^- lO"*"* ergs s'K 

M06 interpreted the result in Sazonov et alJ (l2006l hereafter 
S06) to be 1.0 ±0.3 x 10^^ ergs s'^ for 1032-7-34 grgs s'^ 
using Eq. 5 in S06 and Eq. 7 in M06, and claimed their result 

" M06 assume F = 1.5 for the X-ray spectra of the sources in the 2° X 1° 
region around the GC 
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(5 ±2 X 10^^ ergs s"' Mq) is consistent with S06. However, 
the result in S06 is calculated in the 2-10 keV range and M06 
in 0.5-8 keV. In the 2-8 keV range, the result in M06 be- 
comes 3.3 ± 1.3 X 10^^ ergs s"' Mq under their assumption 
of r = 1 .5 and A^h = 6 x 10^^ cm"^. Similarly, the result in S06 
is scaled to be 9 ± 3 x 10^^ ergs s"' Mq in the same energy 
band. So there is a hint of mismatch in the results between 
M06 and S06 if one follows the scaling in M06, but also note 
Fig. 9 in S06 shows - 6 ± 2 x lO^*" ergs s"' M"' in the 2- 
10 keV band for > 10^^"^^ ergs s~', which is consistent with 
M06 and lower than the scaling done for S06 in M06. This 
conversion also reveals the result in M06 is consistent with 
our result for the Sgr A* field. Using a similar scaling based 
on Eq. 5 in S06, we get 3 ± 1 x lO^'' ergs s"' M;^' for the X- 
ray em issivity reported by iRevnivtsev. VikhUnin & Sazonovl 
mO% hereafter R07b) in the 2-8 keV range for lO^^-^-^^ 
ergs s~', assuming 3% of tht total emission in the same lumi- 
nosity range based on Fig. 6 in R07b. Due to many different 
underlying assumptions in the above estimates (e.g. the spec- 
tral model parameters, stellar models, etc), it is not easy to 
make a fair comparison among the reported results. The large 
uncertainties make these results appear consistent within 2a, 
but our res ults are at t he lower end of these findings. 

M06 and lMuno eTa l. (2008, hereafter M08) have presented 
the X-ray source distribution in a 2° x 1 ° region around the 
GC. In the case of the logA^-log5 distribution, one of the inter- 
esting results in M06 and M08 is a flatter distribution of the 
X-ray sources in the Arches Cluster and the subsequent ex- 
cess of the X-ray sources near the high end of the flux range, 
compared to the Sgr A* field. We also see a similar cross over 
between the unstacked Sgr A* field and the Arches Cluster at 
^ 1.5 X 10~'^ ergs cm~^ s"' (or at 8 x 10"'^ ergs cm"^ s"' 
with the stacked Sgr A field, out of the range in the Fig. 5a). 
We believe this is a simple statistical fluctuation rather than 
a true representative of the population, arising from a small 
number of sources in the narrow FoV, where a few strong 
sources (^^ 2-4 in the Arches Cluster) skew the shape of the 
whole distribution. In fact, the small number statistics is also 
evident in the jumpy shape of the distributions near the high 
end of the flux range. In Fig. 5 of M06, the excess of the X- 
ray sources in the Arches Cluster above 6 x 10"^ ph cm"^ s~' 
is boosted by excluding the overlapping region between the 
Sgr A* field and the Arches Cluster from the sky coverage 
calculation for the X-ray sources of the Arches Cluster. In or- 
der to minimize the effects due to the small number statistics, 
in our analysis, we focus on the solid line section of the distri- 
butions that contain at least 10 or more sources. In addition, 
quantile analysis results in a slightly higher value (^ 10%) 
of the rate-to-flux conversion factor for the G3 sources in 
the Sgr A* field than the same for the Arches Cluster (see 
Fig. 4), which in turn pushes the logA^-log5 distribution of 
the Sgr A* field relatively higher. As a result, we find the 
slopes of the log^V-log^ distributions of the Arches Cluster 
and Sgr A* fields are consistent and the Sgr A* field contains 
more X-ray sources than the Arches Cluster consistently be- 
low ~ 3 X 10"'^ ergs cm"- s"'. 

In the case of the radial distribution, the excess of the X-ray 
sources in the Sgr A* field with respect to the stellar model 
is about 3(7 above the stellar model if one considers both the 
statistical errors and the ^ 20% systematic errors in the flux 
estimates (about IQa above only with the statistical errors). 
M08 find about 2.5ct excess of the X-ray sources at the GC 
compared to the best fit stellar model. 



5.4. Another source population in Sgr C (or Sgr B2)? 

The 1/9 distribution is roughly consistent with the observed 
radial distribution of the X-ray sources in the GB fields within 
2-3(7 except for the Sgr C field. The 1/9 distribution is a 
merely empirical outcome if the observed X-ray source den- 
sity consists of multiple components. On the other hand, the 
apparent 1/9 distribution along with the similar X-ray spectral 
properties may imply that the GB X-ray population belongs to 
a homegenous component. If true, the excess of the Sgr C and 
Sgr B2 fields can be simply viewed as the presence of another 
source population in these fields in addition to the population 
following the 1/9 distribution. The Sgr C field, like the Sgr 
B2 field, contains molecular HII complexes that host massive 
star formation. These molecular clouds are very luminous in 
hard X-rays, in particul ar with the 6.4 keV neutral iron line 
jMurakami et al.1 12001 alibi) . In our analysis, the estimates of 
the r value in the PL and PL-nFe model for the stacked spec- 
tra of the G3 sources in the Sgr B2 and C fields, are relatively 
lower compared to the rest of the fields, suggesting the pos- 
sibility of another source population with a different spectral 
type that could be related to the star formation. 

6. CONCLUSION AND FUTURE WORK 

In the logA^-log5 distribution of the sources in the GB 
fields, the systematic errors arising from certain assumptions 
of spectral type is usually disregarded due to lack of alterna- 
tive approaches. However they often dominate other system- 
atic errors such as the EB, completeness or even statistical er- 
rors. The quantile analysis allows for a simple, robust method 
to assign a proper spectral type for flux calculation. The tech- 
nique is shown to be reliable in the hard band (> 2 keV) and 
insensitive to selection of the spectral model. In the soft band 
(< 2 keV), where the Galactic extinction has a great influence 
in the spectra, the result can vary drastically depending on the 
assumed spectral model class. Therefore, any results covering 
the soft energy range should be taken with caution. 

The logA^-log5 and radial distribution of the GB fields in- 
cluding the three low extinction Windows show the high con- 
centration of the GB X-ray sources near the GC. The GB dis- 
tribution clearly extends out to ^ 1 .4° (LW) from the GC and 
possibly more. The spectral type of the GB X-ray sources 
appears to be largely consistent across the region under the 
power law model with an iron emission line at 6.7 keV. It is 
possible that one type of source constitutes the majority of 
the GB population, and the estimated X-ray density is con- 
sistent with the majority being magnetic CVs (IPs). The 1/9 
relation of the radial distribution of the X-ray sources appears 
to be simply empirical in comparison with the stellar model 
compositions. Since the gravitational influence of the central 
supermassive black hole is only dominant within a few central 
pc (a few arcmin) at most, there is no compelling physics be- 
hind the 1/9 relation being fundamental out to a few degrees 
from the GC. However, the discrepancy between the stellar 
models and the X-ray distribution, and the apparent homo- 
geneity of the X-ray spectral properties of these X-ray sources 
pose other possibilities: the 1/9 relation may not be so em- 
pirical or at least the X-ray source population contains some 
components that are not easily traceable by the visible stellar 
population. 

The radial distributions of the total X-ray and IR surface 
brightness in the fields match within the given uncertain- 
ties, implying the same origin. The radial distribution of the 
IR surface brightness resembles the distribution of the X-ray 
point sources perhaps better than predicted by stellar distribu- 
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tion models, but a further detailed analysis is required because 
of the sensitive dependence of the IR surface brightness on the 
extinction estimates in the high extinction fields. In the case 
of the X-ray surface brightness, the additional care must be 
taken due to dominant contributions of a few bright sources 
on the X-ray flux in star formation fields such as The Arches 
field. 

While multiwavelength observational campaigns provide 
important clues on the GB X-ray population, the true nature 
of the nature of GB X-ray sources may not be completely re- 
solved due to source confusion a nd high obscuration. A deep 
observation 1 Ms) of the LW jRevnivtsev et alll2009l) . de- 
signed to investigate the nature of the Galactic Ridge X-ray 
emission (GRXE) in the field, is very encouraging for stud- 
ies of the nature of X-ray point sources in the GB. Such a 
deep observation allows a direct detection of iron emission 
lines or X-ray variability in many of the GB X-ray sources, 
with which we can identify the nature of individual sources. 
We note only a handful of sources in the 1 Ms data of the 
Sgr A* field had enough statistics for identification through 
such a direct discovery (M03,M04). But the low extinction in 
the Window fields can be a game changer. For instance, we 
have identified an IP in the BW from the 100 ks observation. 

Assume the identifiable source distribution is proportional to Sth^^^, 
where Sn, gets 10 times fainter, and also assume an additional 20-30% in- 
crease in the probability of catching highly variable X-ray sources based on 
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TABLE 1 

X-RAY POINT SOURCES IN THE SELECTED GB FIELDS 



Field 


Obs. ID 


/ 


b 


"Offset 




'GTI 




Source count'' 








n 


n 


n 




(ks) 


Bx 


Sx 


Hx 


Combined 


BW 


3780 


1.06 


-3.83 


3.93 


0.31 


96 


365 


326 


134 


407 


SW 


4547, 5303 


0.25 


-2.15 


2.12 


0.48 


96 


388 


313 


140 


433 


LW 


5934, 6362,6365 


0.10 


-1.43 


1.39 


0.68 


94 


282 


184 


174 


319 


SgrB2 


944 


0.59 


-0.03 


0.65 


81.2 


97 


279 


126 


224 


363 


SgrC 


5892 


-0.57 


-0.02 


0.51 


52.7 


97 


313 


188 


241 


442 


Arches 


4500 


0.12 


-0.02 


0.19 


52.5 


97 


330 


84 


328 


423 


Sgr A* 


3665 


-0.06 


-0.05 




56.5 


88 


401 


92 


400 


508 


(Stacked^) 












698 


2251 


370 


2316 


2876 



(a) The aim point offs et from Sgr A*, (b) The estimates for the integrated neutral hydrogen column density along the line of the sight (in 
10^" cm"") by Schlegel et al.l d 19981) for the location of the aimpoint. This is only for guiding purpose due to the large uncertainty in the Galactic 
plane fields, (c) The good time intervals (GTIs). The total exposure (i.e. before cleaning) is 100 ks each (750 ks for the stacked Sgr A* field), 
(d) The number of the sources with net count > 1 in the broad band (0.3-8 keV) on the ACIS-I CCDs (0, 1, 2, and 3) in the three detection 
bands (Bx'- 0.3-8 keV, Sx'- 0.3-2.5 keV, Hx- 2.5-8 keV) and the combined unique source list, (e) 14 pointings are stacked, and they are Obs. ID 
242, 2951, 2952, 2953, 2954, 2943, 3663, 3392, 3393, 3549, 3665, 4683, 4684 and 5360. 



TABLE 2 

Catalog of X-ray point sources in the Window and four GB fields 



Source 








Posi. 






Net counts'' 








Quantiles 




Unabsorbed Flux' 


Name 


Field 


R.A. 


Dec. 


EiTor^' 


Oftset'' 


Bx 


Sc 


He 


He 


£^50 


Quaitile 


Group 


He 


(CXOPS J) 




(°) 


(°) 


(") 


(') 










(keV) 


Ratio' 




(10"" ergs cm"- s"') 


180230.4-295647 


BW 


270.626934 


-29.946497 


1.29 


10.05 


123.6(14.3) 


71.5 (10.9) 


53.3 (9.9) 


5.4 


1.89 (0.10) 


1.30 (0.17) 


1 


1.31 (0.24) 


18023 1. 2-295528 


BW 


270.630007 


-29.924698 


2.37 


10.12 


38.6(10.9) 


4.6 (6.9) 


34.8 (8.9) 


3.9 


3.56 (0.29) 


1.68 (0.30) 


3 


1.12(0.29) 


180235.9-295323 


BW 


270.649946 


-29.889846 


3.32 


9.87 


23.6 (9.6) 


23.4 (8.1) 


-0.2 (5.6) 


0.0 


1.05 (0.15) 


2.08 (0.48) 


1 


-0.01 (0.13) 


175404.4-294359 


SW 


268.518385 


-29.733089 


2.56 


9.58 


34.1 (11.2) 


25.3 (9.2) 


5.9 (6.6) 


0.9 


1.34(0.18) 


1.63(0.64) 


1 


0.13(0.14) 


175405.3-294717 


SW 


268.522117 


-29.788307 


2.47 


8.04 


22.2 (8.9) 


19.7 (7.7) 


2.9 (5.1) 


0.6 


1.40(0.24) 


1.61 (0.40) 


1 


0.06(0.11) 


175406.7-294239 


SW 


268.527957 


-29.711050 


1.97 


9.99 


42.2 (12.0) 


20.6 (9.4) 


22.5 (8.0) 


2.8 


2.12(0.54) 


1.08 (0.32) 


2 


0.69 (0.24) 


175051.2-293418 


LW 


267.713518 


-29.571797 


1.02 


8.10 


113.8(13.8) 


36.2 (9.1) 


75.8(10.9) 


7.0 


2.63 (0.20) 


1.25 (0.16) 


2 


2.21 (0.32) 


175052.0-293319 


LW 


267.716827 


-29.555400 


2.92 


8.14 


16.4 (9.2) 


5.2 (6.8) 


9.3 (6.6) 


1.4 


0.97 (5.31) 


0.20 (0.29) 


2 


0.27 (0.19) 


175053.3-293207 


LW 


267.722097 


-29.535548 


2.12 


8.29 


25.0(10.3) 


2.8 (7.3) 


22.4 (7.8) 


2.9 


3.46 (0.38) 


1.74 (0.48) 


3 


0.77 (0.27) 



Notes. — This table shows a part of the complete list, which is available in the electronic edition. 

(a) The 95% positional eiTor radius, (b) The offset from the aim point, (c) The net counts based on the aperture photometrv lHong et aU20Q3) . (d) The signal to noise ratio in the He band. The sources with S/N > 3 are included in the logW-log5 
plot in Fig. 5. (e) 3(£'25 —0.3 keV)/(£'75 —0.3 keV). (f) Based on the PL+FE model using quantile analysis. We do not include the flux estimates in the other bands due to their large uncertainty. See the text for the details. The uncertainties for net 
counts and fluxes are statistical enors. 
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TABLE 3 

SPECTRAL MODEL PARAMETERS FOR THE G2 AND G3 SOURCES 







Quartile 




PL from 


PL+Fe (Ht 


:-a) from 




PL (+Fe He-Q) from 




Field 




Ratio" 


Quantile Diagram 


Quantile Diagram^ 




Spechral Model Fit'' 






£^50 




i 


%22 


i 


%22 


r 


%22 


EW 






(keV) 






(xlO^^ cm" 


2-1 
} 


(XlO^^ cm"^) 




{ V 1 Cr^ cxrT- 


-) (keV) 




unabsorbed hard sources 


(G2) 


















BW 


1.97(2) 


1.02(1) 


1.38(03) 


0.31(03) 


1.42(03) 


0.32(05) 


1.36(2) 


0.26(01) 


0.15(9) 


113. 3/167 


SW 


2.10(5) 


1.03(2) 


1.35(07) 


0.37(06) 


1.38(10) 


0.38(08) 


1.22(3) 


0.25(03) 


_g 


48.7/69 


LW 


2.66(4) 


1.09(2) 


1.28(07) 


0.73(09) 


1.35(07) 


0.79(09) 


0.99(2) 


0.38(02) 


_E 


125.9/148 


absorbed hard sources (G3) 


















BW 


3.22(13) 


1.38(9) 


1.66(37) 


2.20(70) 


1.74(37) 


2.30(70) 


1.22(4) 


1.66(16) 




24.8/26 


SW 


3.39(7) 


1.52(4) 


1.77(23) 


2.90(50) 


1.91(23) 


3.10(50) 


1.58(4) 


2.78(19) 




35.4/22 


LW 


3.48(5) 


1.43(3) 


1.21(10) 


1.95(15) 


1.32(10) 


2.10(20) 


1.30(2) 


1.89(06) 


0.17(8) 


93.1/115 


Sgr B2 


4.75(4) 


1.86(2) 


-0.37(14) 


3.40(80) 


0.25(17) 


5.70(90) 


0.50(1) 


6.20(22) 


0.61(7) 


105.2/154 


SgrC 


4.81(3) 


1.90(2) 


-0.26(10) 


4.8(1.0) 


0.46(21) 


7.8(1.2) 


-0.10(1) 


3.95(18) 


0.38(5) 


172.3/189 


Arches 


4.54(2) 


1.83(1) 


0.14(07) 


4.00(50) 


0.67(14) 


5.70(70) 


0.85(1) 


5.17(12) 


0.66(5) 


319.4/363 


Sgr A* 


4.69(2) 


1.91(1) 


0.31(14) 


6.40(60) 


0.94(14) 


9.00(80) 


1.02(1) 


6.95(12) 


0.46(4) 


324.3/364 


(a) 3(£25 ■ 


-0.3 keV)/(£75- 


-0.3 keV). (b) The parameter estimates based 


on quantile analysis for 


a power law model, (c) The same 


as (b) but with 


a fixed 0.4 keV EW at 6.7 keV (d) 



The parameter estimates by the spectral model fit. (e) the EW of 6.7 keV line, (f) Degrees of Freedom, (g) Due to poor statistics, the spectral fit is done with a power law model without 
an iron line. See g] 
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TABLE 4 

X-RAY SOURCE AND STELLAR DENSITY" 





X-ray 


Source^ 






Stellar Model A'' 






Stellar Model B" 




Surface 


Volume 


Star volume 


Star surface 


X-ray to 


X-ray 


Required CV 


Star volume 


Field 


density 


density 


density 


mass density 


stellar mass 


to stars 


to stars'^ 


density 




(deg--) 


(IQ-^ pc-') 


(pc"') 


(IO^Mq deg--) 


(lO-'M^) 


(10-') 




(pc-^) 


Field to field Compari: 


son 
















BW 


25(45) 


3.2(5.7) 


2.1 


0.5 


5.1(9.1) 


1.6(2.8) 


2.6(4.6) 


3.2 


SW 


39(47) 


4.9(6.0) 


5.3 


1.3 


3.1(3.7) 


0.9(1.1) 


1.6(1.9) 


6.7 


LW 


203(64) 


25(8.1) 


7.6 


1.8 


11(3.5) 


3.4(1.1) 


5.6(1.8) 


8.3 


Sgi- B2 


582(93) 


73(11) 


41 


10 


5.7(0.9) 


1.8(0.3) 


2.9(0.5) 


8.8 


Sgi-C 


902(113) 


113(14) 


45 


11 


8.1(1.0) 


2.5(0.3) 


4.1(0.5) 


*41 


Arches 


1624(153) 


204(19) 


52 


13 


12(1.2) 


3.9(0.4) 


6.5(0.6) 




Stacked Sgr A* 


3100(187) 


389(23) 


68 


17 


18(1.1) 


5.7(0.3) 


9.5(0.6) 




100 ksSgiA* 


2117(176) 


266(22) 


68 


17 


12(1.1) 


3.9(0.3) 


6.5(0.5) 




By the fit to the radial distribution' 
















with stacked Sgr A* (Fig. 5b or 6a) 










10.3(0.5) 


3.2(0.2) 


5.3(0.3) 


with lOOks Sgr A* (Fig. 6b) 










9.2(0.5) 


2.8(0.2) 


4.7(0.3) 


By the fit with the freed relative normahzation of p^^ 














with stacked Sgr A* (Fig. 6c) 










7.1(0.6) 


2.2(0.2) 


3.6(0.3) 


with lOOks Sgr A* (Fi 


g. 6d) 










7.5(0.6) 


2.3(0.2) 


3.8(0.3) 



(a) Assumes the hard X-ray sources {S > 1.5 X 10 ergs cm ^ s ' in He) or stars seen within the 17' X 17' FoV are mainly (^ 80%) from 6-10 kpc distance, (b) Using the PL+Fe 
model. We subtract the expected AGN numbers, 145 deg~^ from the Window fields and 107 deg~^ from the high extinction fields, (c) The composite model in M06 and references 
therei n. The model gives the stellar mass density in the unit of Mq pc~^ , and we assume the local value of 0. 144 stars pc~^ and 0. 04 pc~^ to get th e star number density i Robin et al] 
1 20031). Th is relation should be good for the bulge in the case of CVs and active binaries, but perhaps not good for young stars iSazonov et aU20Q^ . (d) The stellar density model by 
iPicaud & Robin i 2004) for the Galactic disk and outer Galactic bulge. The model does not include a central nucleus, so the (italic) values for the Sgr B2, Sgr C, Arches, and Sgr A* 
fields are not reliable. See j |5.1| (e) The required CV to star density to explain the hard GB X-ray sources by IPs. We assume that IPs are 5% of all CVs (e.g. ~ 2-8% in R06) and that 
about 12% of them are detected above 10^" ergs s~^ (e.g. ~ 10-16% in R06), which conesponds ~ 1.5 X 10"'"^ ergs cm~" s~^ for the sources near the GC (see text), (f) The relative 
normalizations among stellar model components are fixed as given in Section 4.3. (g) The relative normalization of the central spherical cluster component of the stellai" distribution is 
allowed as a free parameter and fitted as well. 



TABLE 5 

Total X-ray and IR surface brightness (see also Fig.[7]i 







X-ray" 






IR 




Field 


2-8 keV (the He band) 




3.6 jim 






(10"'" erj 


;s cm " 


deg-2) 




(MJy/sr) 






PL 




PL+Fe 


observed 


unabsorbed (PL)" 


unabsorbed (PL+Fe)^ 


BW 


0.34(0.02) 




0.35(0.02) 


7.4(0.4) 


8.2(0.4) 


8.2(0.4) 


SW 


0.48(0.01) 




0.49(0.01) 


15(1) 


17(2) 


17(2) 


LW 


LlO(O.Ol) 




1.13(0.01) 


18(1) 


24(1) 


24(1) 


SgrB2 


5.09(0.05) 




5.42(0.05) 


26(6) 


74(18) 


150(36) 


SgrC 


14.9(0.1) 




16.6(0.1) 


35(3) 


152(13) 


384(33) 


Arches 


30.9(0.1) 




33.4(0.1) 


58(12) 


199(40) 


337(69) 


Sgr A* 


23.93(0.02) 




27.72(0.02) 


68(16) 


493(110) 


1100(250) 



(a) Assume the spectral model using the G2 group for the Windows field and the G3 group for the rest in Table [3| The error is of the 
statistical origin. Note the expected cosmic X-ray background ( CXB) is (1.6- 1.7 ± 0.2 ) x 10"" ergs cm"^ s-' deg-' for the Windows fields 
and ( 1 .3 - 1 .4 ± 0.2) X 10-" ergs cm-^ s-' deg-^ for the rest i Hickox & Markevitch II2006I) . See the fcl (b) The Na estimate by PL (line B in 
Fig.|7J. (c) The A'h estimate by PL+Fe (line A in Fig.|7j. 



TABLE 6 

Specific Luminosity of X-ray point sources 







Energy 


Luminosity 


Scaled for 




Source 


Reported 


Range 


Range 


10" '--" ergs S-' in2-8keV 


Studied Fields 




(10-'' ergs S-' Af^) 


(keV) 


(ergs s-') 


(lO-"* ergs S-' M^) 




This Study 


0.5-2.8 


2-8 




0.5-2.8 


7 GB fields (100 or 750 ks) 


M06 


5.0 ± 2 


0.5-8 


J032.7-34 


3.3 ± 1.3 


2° X 1° around the GC (100 or 2 X 12 ks ) 


S06 


45 ±9 


2-10 


J027-36 


9±3' 


the local solar neighborhood 










6±2'' 




R07b 


77 ± 39 


2-10 


JQ30.3-32.3 


3± I' 


the Sgr A* field (1 Ms) 



M06 - iMuno et al. ( 2606). S06 - ISazonov et alj|2006h . R07bTRevnivtsev. Vikhlini n & Sazono\< 12001 
(a) by the scahng by M06. (b) Fig. 9 in S06. (c) assuming 3% of the total emission ((4 ± 2) X lO" ergs s-' ) is from > 10^^ ' ergs s-' , based on Fig. 6 in R07b. 
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0.01 0.10 1.00 10.00 100.00 0.01 0.10 1.00 10.00 100.00 0.01 0.10 1.00 10.00 100.00 
d, (offset/quod, sum of pos. error) (offset/quod, sum of pos. error) d^ (offset/quod, sum of pos. error) 



Fig. 1 . — Cross correlating X-ray sources detected in three different energy bands: tlie number of pairs as a function of tlie relative separation (dr, see text) of 
potentially identical sources (closest pairs) among the three detection bands. The plots show the results from the 100 ks observation of the Stanek Window (a), 
the Sgr A* field (b) and 750 ks of the Sgr A* field (c). The bimodal shape is due to the mixture of the true and random matches in the distribution. The (blue) 
dashed-dotted lines show the results after introducing an arbitary global offset (1' in both R.A. and Dec.) among the three band detections, which illustrate the 
distributions of the random matches. 
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Fig. 2. — Quantile diagrams (0.3-8 keV) of the X-ray sources with S/N > 3 in the GB fields overlaid with grids for a simple power law model (PL, solid lines, 
power law index T = 0, 1, 2, 3 & 4, %22 = 0-01, 0.1, 0.4, 1 & 10), a power law plus an iron line model (PL-l-Fe, dashed, at 6.7 keV with 0.4 keV EW), and 
thermal Bremsstrahlung model (TB, dotted, kT = 0.2, 0.4, 1, 2, 4 & 10 keV, %22 = "-Ol' O-L 0-4, 1 & 10, only shown in the bottom-right plot for clarity). The 
energy quantile Ey corresponds to the energy below which j(% of the counts are detected. The (red) crosses are for the relatively soft sources (S/N > 3 in Sc, but 
not in He), the (blue) 'x's for the hard sources (S/N > 3 in He, but not in Sc), the (black) triangles for the bright sources (S/N > 3 in both 5c and He), and the 
(orange) dots for the faint sources (S/N > 3 only in Be). The (grey) lines from (-0.5,1.3) divide each diagram into the soft (Gl), medium (G2) and hard groups 
(G3). 0s mark the quantiles of stacked photons in each group. 
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Fig. 3. — The stacked spectra of the G3 sources (net counts < 1000) in LW with the PL+Fe XXV He-a (6.7 keV line) fit. The estimated EW of the line is 
0.17±0.08keV. 
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Fig. 4. — Comparison of the rate-to-flux conversion factor for the three quantile groups of sources in the three energy ranges. The conversion factor in the He- 
band (2-8 keV) is robust ( < 20-30% variation), in the Sc band (0.5-2 keV) it is very unreliable (up to more than a factor of 10), and in the medium energy range 
(1-7 keV), there are significant variations (up to ~ 100%). 
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Fig. 5. — The \ogN-\ogS (a) and the radial (b) distributions of the X-ray sources in the GB fields. Fluxes are computed under the PL+Fe model (6.7 keV line 
with 0.4 keV EW) and the distributions include the sources with 5/^ > 3 in the He (2-8 keV) band. The numbers in the legend of (a) are the slopes (a) and their 
error of the power law fits (N oc S^" ) to the solid sec tion of the lo°A^-log5 distributions (see the text for the definition of the solid section). The (orange) solid line 
is the active galactic nuclei (AGN) distribution from lKim et in (2007) seen in the low extinction fields and the (orange) dashed line is the same corrected for the 
extinction to the GC (A'jj =6x 10^^ cm"^, see the text). The radial distribution shows the number density of the X-ray sources with S > 1.5 X 10"''' ergs cm"^ s"' 
(marked by the vertical strip in the left panel) under two spectral models (solid black for PL-l-Fe and dotted red for PL), compared with the stellar distribution 
(solid green) and the 1/9 distribution (blue dashed). The.v-axis in the radial distribution is the average offset value of the sources in each field. 
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Fig. 6. — Comparison of the X-ray and stellar distributions: the best model fit of the stellar to X-ray distrubtion with (a) the stacked Sgr A* data, and (b) the 
100 ks Sgr A* data, and (c & d) the same model fit but with the freed normalization parameter of the central spherical cluster component (pa)- The excess of the 
central spherical cluster component (p^i) needed for the best model fit is 8.4 or 4.5 X higher than the original model, depending on the Sgr A* data set. 
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Fig. 7. — Comparison of the X-ray and the IR surface brightness. The X-ray surface brightness (black solid line) is the total unabsorbed Hq flux minus the 
expected CXB flux in the Chandra ACIS-I FoV of each field assuming the quantile average spectral model (PL+Fe) of the G2 (for the Windows fields) or G3 group 
(for the rest). The infrared surface brightness (blue solid) is the average flux in the Chandra ACIS-I FoV from the 3.6 /.tm band images by Spi'fter/GLIMPSE. To 
calculate the unabsorded IR flux consistently with the X-ray flux, we assumed the estimates by the quantile analysis in Table[3|(G2 for the Windows fields 
and G3 for the rest). We use the N^i estimates by both models (A for PL and B for PL+Fe). We also assumed A3 Sfim = 0.5 Aks I N ishiv ama et al. 2009). The 
errors of the X-ray brightness are the quadratic sum of the statistical and systematic errors and the errors of the CXI3 estimates (see j|5.21 . The errors of the IR 
brightness are based on the statistical fluctuation of the brightness in the four ACIS-I fields, but the dorminant errors in the high extinction fields come from the 
uncertainty in the extinction estimates as seen in the large difference of lines A and B. For easy comparison, we overlay the distributions of the X-ray source 
density (black dotted) and what the stellar model predicts (green dotted) without the AGN component using an arbitary normalization factor. 



